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ABSTRACT 

We present an X-ray spectral analysis of 126 galaxies of the 12 micron galaxy sam- 
ple (12MGS). By studying this sample at X-ray wavelengths, we aim to determine the 
intrinsic power, continuum shape and obscuration level in these sources. We improve 
upon previous works by the use of superior data in the form of higher signal to noise 
spectra, finer spectral resolution and a broader band pass from XMM-Newton. We 
pay particular attention to Compton thick AGN with the help of new spectral fitting 
models that we have produced, which are based on Monte-Carlo simulations of X-ray 
radiative transfer, using both a spherical and torus geometry, and taking into account 
Compton scattering and iron fluorescence. We use this data to show that with a torus 
geometry, unobscured sight lines can achieve a maximum equivalent width (EW ) of 



the Fe Ka line of 150 eV, originally shown by Ghisellini, Haardt fc Matt| ( |1994[ ). In 
order for this to be exceeded, the line of sight must be obscured with A'h> 10^"^ cm~^, 
as we show for one case, NGC 3690. We also calculate flux suppression factors from 
the simulated data, the main conclusion from which is that for A^h^^ 10^^ cm~^, the 
X-ray flux is suppressed by a factor of at least 10 in all X-ray bands and at all red- 
shifts, revealing the biases present against these extremely heavily obscured systems 
inherent in all X- ray surveys. Furthermore, we confirm previous results from [Murphy] 



&: Yaqoob (2009) that show that the reflection fraction determined from slab geome- 
tries is underestimated with respect to toroidal geometries. For the 12 micron selected 
galaxies, we investigate the distribution of X-ray power-law indices, finding that the 
mean (< F >= 1.90to;o7 ^'^'^ = 0.3llg;Q5) is consistent with previous works, and 
that the distribution of F for obscured and unobscured sources is consistent with the 
source populations being the same, in general support of unification schemes. We de- 
termine a Compton thick fraction for the X-ray AGN in our sample to be 18 ± 5% 
which is higher than the hard X-ray (> 10 keV) selected samples. Finally we find that 
the obscured fraction for our sample is a strong function of X-ray luminosity, peaking 
at a luminosity of ~ iQ'*^"^'^ ergs~^. 
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1 INTRODUCTION 

Studying AGN at X-ray wavelengths allows an assessment 
of the intrinsic source properties, such as the measurement 
of the continuum and the source power, due to the ability of 
X-rays to penetrate and measure large columns of obscuring 
material. This makes X-ray observations essential for test- 
ing AGN unification schemes ( [Awaki et al.^^1991) which in- 
voke obscuration to explain the perceived difference between 
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Seyfert Is and Seyfert 2s ( |Antonucci|1993[ |Urry fc Padovani| 
1995 1. The Ah and T distributions for AGN are also an im- 



portant ingredients for synthesis models of the cosmic X-ray 
background (XRB, |Gilh, Comastri fc Hasinger||2007[ ). 

Most of our knowledge of the X-ray properties of lo- 
cal AGN comes from X-ray selected samples, such as the 
[Piccinotti et al.| (|1982[) sample de rived from HE AO 1 ob- 
servations (e.g. Turner fc Pounds 1989 Nandra fc Pounds 
19941. X-ray selection of AGN is in principle very effective 



as in this band the nucleus predominantly outshines the 
host galaxy, and can penetrate obscuration. The hard X- 
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ray telescopes aboard INTEGRAL and Swift and their large 
sky coverage have allowed for the compilation of hard X- 
ray selected AGN catalogues. 'Beckmann et al. (20061 and 
[Beckmann et al. (2009J have presented catalogues for INTE- 
GRAL selected AGN, which total 199 after 5 years of obser- 



vations, and Tueller et al. (20081 and Tueller et al. (20101 



have presented 229 Swift/BAT selected Seyfert galaxies af- 
ter 22 months of observations. However, the discovery from 
these surveys seems to be that the fraction of Compton thick 
AGN detected is smaller than expected from XRB studies 



(Gilli, Comastri & Hasinger 20071. This points to the em- 



pirical fact that even the X-ray band is biased against the 
most obscured, Compton thick, AGN. X-ray selection at en- 
ergies greater than 10 keV is less affected by Compton thick 
obscuration than 2-10 keV selection though. 

An alternative to X-ray AGN selection is mid-infrared 
(MIR) selection where the primary radiation of the AGN is 
re-emitted after having been reprocessed by hot dust. The 
extended IRAS 12 micron galaxy sample (12MGS - |Rush,| 



Malkan fc Spinoglio|1993l RMS) is a sample of 893 MIR se- 
lected local galaxies which contains a high fraction of AGN 
(13% at the time of publication, RMS). The sample is taken 
from the IRAS Faint Source Catalogue, version 2 (FSG-2) 
and imposes a flux limit of 0.22 Jy, including only sources 
with a rising flux density from 12 to 60 microns (to exclude 
stars) and with a galactic latitude of Js25 deg. Being 
selected in the mid-infrared (MIR) this sample is also rela- 
tively unbiased against absorption, low luminosity systems 



and 'hidden AGN'. Spinoglio & Malkan (19891 showed that 



a wide variety of AGN types emit a constant fraction of 
their bolometric luminosities at this wavelength, and fur- 
thermore shown to be true for star forming galaxies as well 



by [Spinoglio et al.| ( |1995| ). The 12MGS should be therefore 
also representative of the true number of different active 
galaxy types in relation to each other. 

It is possible, though, that even at 12 |j.m, the AGN 
emission may be suppressed in the most heavily obscured 
nuclei. Pier fc Krolik] ([l992) present radiative transfer mod- 
elling for the infrared emission from dust tori for varying 
torus orientations. The most significant result of this study 
was that throughout the infrared, the emission in the edge- 
on direction is weaker with respect to the face-on direction, 
and at 12 ^.m, the difference can be up to several orders 
of magnitude for optically thick tori. These models suggest 
that AGN selection at 12 microns is biased against heavily 
obscured nuclei. However, recently this was tested observa- 



tionally by Horst et al. (20081. They used high resolution 



12.3 nm imaging of a sample of local Seyfert nuclei, with 
X-ray data to explore the ratio of the MIR to intrinsic X- 
ray nuclear luminosities in order to detect any difference 
between the obscured nuclei (Seyfert 2s) and unobscured 
nuclei (Seyfert Is). As the X-ray luminosities have been 
corrected for absorption, Lx should represent the intrin- 
sic emission from the nucleus. For smooth tori modelled by 
[Pier &: Krolik| ( [l992| ), obscured systems should then present 
lower Lmir/I/x ratios than unobscured systems due to the 
suppressed MIR flux seen in these systems. A difference of 
an order of magnitude for a difference in A''h of between 
~ 10^'' and ~ 10^'' cm is predicted by the torus mod- 
els, and should be measurable. However, this is not seen 
observationally, as |Horst et al.| ( [2008[ ) find no significant dif- 
ference between Seyfert Is and Seyfert 2s, or even Compton 



thick nuclei. They argue from this that the torus may not 
be smooth, and instead clumpy, where the torus is seen to 
be optically thin, but the clumps are optically thick, though 
with a small volume filling factor. Most importantly for this 
work though, it suggests that 12 micron selection is in fact 
not heavily biased against obscuration. This therefore makes 
the 12MGS the ideal parent sample for study of obscuration 
and unification using X-ray observations and is as such the 
subject of this paper. 

The 12MGS has previously been studied at X-ray wave- 



lengths by Barcons et al. ( 1995 1 who used hard X-ray obser 



vations by HEAO 1, and at soft X-rays by Rush et al. ( 1996 \ 
with ROSAT data. [Barcons et al.| ( |1995[ ) use the X-ray and 
IR data sets to show that most of the local X-ray emissivity 
originates from Seyfert galaxies and that most of the local 
2-10 keV X-ray luminosity function between 10*^ and 10*^ 
ergs~^ can be accounted for by 12 micron emitting AGN. 
They also test the unification scheme and come to the con- 
clusion that it must be modified to include the reduction of 
the covering fraction of the torus with increasing source lu- 



minosity, an observation first suggested by Lawrence ( 1991 \ 



which has come to be known as the 'receding torus model'. 



Rush et al. (1996) conducted a spectral analysis of the 0.1- 



2.4 keV ROSAT data.. For spectra with enough counts, they 
fit absorbed power-laws, finding a median soft X-ray spectral 
index of 2.3 for all Seyfert types. 

In this paper, we present an X-ray spectral analysis of 
the galaxies of the 12MGS with an observation by XMM- 
Newton, which consists of a heterogeneous mix of Seyferts, 
LINERs and star- forming galaxies. We focus on the X-ray 
spectral properties of the sources without regard for their 
optical spectroscopic types and use X-ray luminosity to se- 
lect unambiguous AGN. In a companion paper though, we 
continue our investigation into this sample combining the 
X-ray data from this paper with optical spectroscopic and 
infrared data. In Section [2] we present work on Monte-Carlo 
simulations of X-ray radiative transport which we use in 
our spectral analysis. Section[3]describes our spectral fitting 
methods and Section |4] gives our spectral fitting results. Fi- 
nally we discuss all of these results in Section |5] and present 
our conclusions in Section [S] 



2 NEW X-RAY SPECTRAL MODELS FOR 
HEAVILY OBSCURED SOURCES 

In our X-ray spectral analysis of the 12MGS, we pay partic- 
ular attention to whether sources are Compton thick, as the 
Compton thick fraction is particularly crucial for the XRB 
models and knowledge of the accretion history of the uni- 
verse. Compton thick AGN are most commonly identified 
in the 2-10 keV X-ray band by the measurement of a flat 
spectrum, where F < 1.0, and the equivalent width (EW) of 
the iron Ka line has been measured to be high (EW>1 keV, 
19961. This is typical of a reflection spectrum. 



Matt et al. 



produced by X-ray scattering from an optically thick sur- 
face, such as the accretion disk or wall of the torus in AGN. 
It is indicative of heavy obscuration of the central source 
as the reflection is dominating over the transmission of the 
intrinsic power-law through the obscuring material. 

In some cases, the transmission of the intrinsic power- 
law can be detected in the 2-10 keV band for these heavily 
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obscured sources. However, due to the extreme suppression 
of flux in this band at these high column densities, this is 
very difiicult. Also at the high optical depths in these Comp- 
ton thick sources the calculations for modelling the trans- 
mission spectrum become non-linear due to multiple scat- 
terings. Thus, models describing simple attenuation of the 
spectrum by absorption and scattering become invalid. 

Therefore, accurate models that correctly account for 
Compton scattering, especially multiple scatterings, are nec- 
essary when fitting the spectra of these sources. Currently, 
the best established model within xspec for correctly de- 
scribing X-ray reprocessing by photoelectric absorption and 
Compton scattering is plcabs. This model describes the X- 
ray transmission from an isotropic source at the centre of a 
uniform spherical distribution of matter. It is valid up to en- 
ergies between 10 and 18.5 keV and up to column densities 



of 5x10 cm (Yaqoob 19971. Unfortunately plcabs has a 



limited range of validity and does not include line emission, 
modelling the continuum only. Recently, two works by |Ikeda] 
et al.| (|2TO9| and [Murphy fc Yaqoob| (|2009[) have produced 



improved Monte-Carlo models based on a toroidal obscur- 



ing medium similar to calculations by Ghisellini et al. ( 1994 1 
which include iron line emission. 

Here we have carried out our own Monte-Carlo simu- 
lations based mainly on the methods of |Nandra fc George| 
( 1994 NG94) , but with a point-source of X-rays located at 



the centre of both a spherical and toroidal distribution of 
gas emitting a power-law spectrum. These simulations fol- 
low the propagation of X-ray photons which interact with 
the medium via Compton scattering or photo-electric ab- 
sorption, which can then result in fiuorescent emission. The 
type of interaction is determined by the energy dependent 
interaction cross sections. We calculate the Compton scat- 
tering cross section described in Rybicki & Lightman ( 1986 1 



and calculate the photo-electric absorption cross sections us- 



ing data from Verner et al. (19961, using abundances from 



Anders & Grevesse ( 1989 1 and fluorescent yields from Bam- 



bynek et al. (19721. The medium is assumed to be of con- 



stant density and there is no dependence of the line of sight 
column density on the inclination angle of the observer, with 
the exception of the unobscured sight lines in the torus ge- 
ometry which have a zero column density. Fig. [l] illustrates 
the geometry used in the torus simulations, where the shad- 
ing only serves to differentiate the different opening angles 
used, and does not correspond to density. 

The models importantly also include line emission, most 
significantly iron Ka (6.4 keV), but also iron K/3 (7.06 keV) 
and Ka emission from several other elements described in 
Table [l] As mentioned before, the EW of the iron line is 
strongly dependent on absorption, as is the shape of the 
line's Compton shoulder. The Compton shoulder feature is 
seen on the red side of the line, and is a result of Comp- 
ton down scattering of the first order line and depends on 
the scattering optical depth. In these models, the flux of the 
line and its spectral shape are inextricably linked to A'h and 
therefore it provides additional constraints on this parame- 



ter when spectral fitting. Yaqoob & Murphy ( 2010 1 present 
a detailed examination of the Compton shoulder's proper- 
ties, showing that it is also dependent on spectral shape 
and geometry of the absorber, and hence contains valuable 
information given that it is well resolved. 

We have made several improvements on the previous 




Figure 1. Geometry of the biconical torus distribution of mat- 
ter used. The shading only serves to differentiate the different 
opening angles used, and does not correspond to density. 



work by NG94, many with the aim to improve computa- 
tional times and spectral quality. These include the use 
of a technique described by Lucy (20021, in which radia- 



tion in Monte-Carlo simulations is represented by indivisible 
monochromatic 'packets' of energy which propagate through 
the system in the same way as the photons were simulated 
previously. This technique allowed us to inject any number 
of quanta at each energy without sacrificing the power-law 
relation. Each packet is assigned a 'weight' dependent on the 
width of the energy bin from which is it injected, the number 
of packets being injected at that energy and the energy-flux 
relation taken from the power-law. 



weight = 



AE X F{E) 



A^ packets 

The weights of the emergent packets are then binned 
at the end in the same way. This technique has the ad- 
vantage that it is possible to vary the number of quanta 
released depending on how detailed the spectrum needs to 
be at that energy. It is especially useful for increasing the 
signal to noise of the output spectrum around the iron K- 
shell regime (6-8 keV). Valuable computational time may 
also be saved by decreasing the signal to noise in the fea- 
tureless parts of the spectrum, above 10 keV for example. 
We also made an improvement to the original method by 
introducing variable elemental abundances (with respect to 
the abundance of Hydrogen). The abundance of iron is a 
separate parameter also, as the modelling of the iron lines 
are particularly sensitive to this. 

We have constructed table models for use within xspec, 
which we call trans for the spherical distribution of matter, 
and torus for the toroidal distribution of matter. The pa- 
rameters of the models are the line of sight column density 
of the neutral material (10^° cm"^ ^ Ah ^ 10^*^ cm~^) and 
the photon index of the intrinsic power-law (1 ^ F ^ 3). For 
trans the abundance of iron (0.1 solar ^ Feabund ^ 10 so- 
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Element 


Abundance 


Line energy 


Edge energy 


Fluorescent 




w.r.t Hydrogen 


(keV) 


(keV) 


yield, Y 




(xlO""') 








Carbon 


3.63 


0.28 


0.29 


0.0025 


Oxygen 


8.51 


0.53 


0.54 


0.0086 


Neon 


1.23 


0.85 


0.87 


0.0183 


Magnesium 


0.38 


1.25 


1.31 


0.0303 


Silicon 


0.355 


1.74 


1.85 


0.043 


Argon 


0.0363 


2.96 


3.20 


0.122 


Calcium 


0.0229 


3.69 


4.04 


0.168 


Chromium 


0.00484 


5.41 


6.00 


0.283 


Iron 


0.468 


6.40, 7.06 


7.12 


0.342 


Nickel 


0.0178 


7.48 


8.35 


0.432 



Table 1. Elemental data used in the Monte-Carlo simulations. Abundances are from [Anders &: Grevesse] | |1989[ l, and fluorescent yields 
are from [Bambynek et ah] | |1972[ ) 



lar) and the total elemental abundance (0.1 solar ^ abund. 
^ 10 solar) are included as parameters, whereas for torus, 
the opening angle of the torus (25.8° ^ 9tor ^ 84.3°) and 
the inclination angle of the torus (18.2° ^ ^ 87.1°) are 
included. Both models have validity in an energy range of 
0.1 keV to 320 keV. 



2.0.1 Fe Ka Line Equivalent Width Calculations 

We also use the results from our Monte-Carlo simulations to 
make predictions of the equivalent width of the iron Ka line 
for both geometries given a set of model parameters. In the 
simulations the photon packets produced by a simulated iron 
K-shell fluorescence are flagged and summed as they exit 
the simulation. In this way, we can calculate the equivalent 
width of the line produced, by the sum of iron Ka photon 
flux divided by the continuum level at 6.4 keV and the width 
of the energy bin at this energy. The calculation represents 
the EW of the entire line, including the Compton shoulder. 

Figure [2] shows how the EW varies with A^h, E and 
the iron and elemental abundances for the spherical geome- 
try. The EW of the line increases linearly versus A''h up to 
T ~ 0.1, where the relationship increases faster than linearly 
beyond that. This is due to the suppression of the continuum 
at 6.4 keV at high column densities, against which the EW 
is measured, but possibly also due to the greater number 
of photons at and above the 7.1 keV iron K shell threshold 
having been down-scattered from higher energies. Increasing 
the iron abundance will as expected increase the EW of the 
line due to the greater number of iron atoms per hydrogen 
atom to produce the transition. The abundance of the other 
elements also effects the EW of the iron line, though not as 
significantly. The EW is seen to increase from low elemental 
abundance to ~ 0.3 x solar abundance. This is also expected 
as an increase in the elemental abundances will cause some 
continuum suppression at 6.4 keV, leading to an increase in 
the EW of the iron line. This trend flattens off for iVH< 10^^ 
as the elemental abundance increases beyond solar though, 
but for higher TVh values, the EW decreases above ~ 0.3 x 
solar abundance. Again though, at these extreme A^h values 
the continuum against which the EW is measured is very 
suppressed so it is difficult to know if this behaviour is real. 

Figure[3]shows how the EW of the iron line changes with 



A^H in the torus model for different torus parameters. The 
left plot shows this relationship for equatorial sight-lines, 
where the different lines show different opening angles. The 
smaller the opening angle, and hence the larger solid angle of 
the torus, the larger the EW of the iron line. The EW reaches 
a maximum here of 4 keV, 2.1 keV and 1.6 keV for torus 
opening angles of 9tor~30° , 60° and 75°. The middle panel 
shows the same relationship, but for the polar sight-line. 
The EWs here are systematically lower due to the greater, 
unabsorbed continuum. The EWs reach a maximum of 150 
eV, 96 eV and 63 eV for 30°, 60°and 75°opening angles. 
Due to this dependance on column density and the opening 
angle of the torus, it may be possible to make inferences 
about the circum-nuclear matter distribution in AGN from 
the EW of the iron line, even for unobscured sources. For ex- 
ample, unobscured sources with an EW greater than ~ 150 
eV must have a very narrow opening in the torus, or else 
super-solar iron abundance. Alternatively, the source may 
not be unobscured at all. It may be very heavily obscured 
with the intrinsic emission extremely suppressed below 10 
keV, and only a weak scattered or galactic component visi- 
ble. We use these data later when determining the Compton 
thick nature of our sources. Finally the right panel in this 
figure shows the variation of EW for a fixed torus opening 
angle of 60°, and varying viewing angles, and shows that for 
unobscured sources, the viewing angle does not affect the 
EW of the iron line greatly. 



2.0.2 X-ray flux suppression 

We also use the results from our new models to calculate 
expected fiux suppression factors for different X-ray bands, 
as a function of absorption and for different redshifts. This 
information is useful when estimating the A^h, given only an 
observed fiux and a predicted intrinsic flux, or when con- 
sidering the biases against heavily obscured sources in X- 
ray selected AGN samples. We do this for both spherical 
and toroidal distributions of matter with solar abundances, 
r = 2 intrinsic power-laws and for the case of the torus, a 
60° opening angle and 90° viewing angle. Figure |4] gives the 
fiux suppression factor at different column densities and in 
different bands. The fiux suppression factor is calculated as 



XMM survey of the 12MGS 5 




10^^ 10^^ ir 0.1-1.0 0.1 1.0 10.0 

N[j (cm ) Ye atund/solar elemental abund/solar 

Figure 2. Predictions of the equivalent width of the iron Kq line for the spherical distribution and varying parame- 
ters. Left panel - EW vs Nn with the three lines representing F = 1,2 and 3 in ascending order for solar abundances. 
Middle panel - EW vs iron abundance (in units of solar abundance) with different lines for varying A^h labelled 
in the figure as logio(Arjj/cm~^) and fixed r=2. Right panel EW vs total elemental abundance (in units of solar 
abundance) with different lines for varying N-^ labelled in the figure as logio(.'VH/cm~^) and fixed r=2. 
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Figure 3. Predictions of the equivalent width of the iron K« line for the torus distribution with varying torus 
opening angles and viewing angles. Left panel - EW vs A''h for equatorial viewing angles (78-90°), and varying torus 
opening angles, 9tor- Middle panel - EW vs A^h for polar viewing angles (0-37°), and varying torus opening angles, 
dtor\ Right panel - EW vs ATh for different viewing angles, di, with a fixed torus opening angle, 5toT-=60°; All with 
fixed r=2. 



the intrinsic, unabsorbed energy flux in a particular band 
divided by the flux observed in that band. 

The attenuation of the intrinsic emission by the absorb- 
ing material is evident as the Nyi increases. For the spherical 
case, at a redshift of zero and for energy bands above 10 keV, 
the effect of this attenuation is negligible (< 2) for A^h up 
to ~ 4 X lO^* cm~^, however, at A^H^lO^^'cm"^, the flux in 
these bands is suppressed by a factor of 10. This relation- 
ship is roughly independent of redshift for these bands. For 
the 2-10 keV band, locally, the flux is suppressed by a fac- 
tor of 10 at Nyi=& X 10^"^ cm~^. As redshift increases, the 
ATh at which the flux suppression factor reaches 10 increases 
due to the photo-electric absorption cut-ofF being redshifted 



through the band. The A^h at which the flux suppression 

factor is 10 for rcdshifts one, two and three is 3 x 10^* , 
7 X lO^"* and lO^'' cm~^ respectively for the 2-10 keV band. 

Flux suppression factors calculated from the torus ge- 
ometry behave similarly to the spherical case until ~ 4 x 10^^ 
cm~^ where they flatten off, rather than continuing to in- 
crease as they do for the spherical geometry. This is due 
to reflected X-ray flux from the torus which is present re- 
gardless of the level of obscuration seen directly through the 
torus. The flux suppression factor flattens off at ^ 100 for 
the 2-10 keV band at z=0, however, this improves to ~ 30, 
~ 15 and ~ 12 for z=l,2 and 3 due to this band observ- 
ing higher energy rest-frame X-rays which are less effected 
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by photo-electric absorption that dominate in the 2-10 keV 
band at rest-frame energies. For the > 10 keV bands, the 
flux suppression factor levels off at ~ 10 at all redshifts. 
At z=3, flux suppression in the 2-10 keV band is less than 
that seen in the 20-100 keV (INTEGRAL) and 14-195 keV 
(Swift/BAT) bands. This is most likely due to the Compton 
hump, produced by the 'down-scattering' of higher energy 
X-rays and peaks at rest-frame ~ 30 keV. At z=3 it appears 
directly in the 2-10 keV band, whereas the higher energy 
bands only view the tail of the hump. 

For both geometries, all redshifts and X-ray bands, flux 
suppression is > 10 at A^H=10'^^cm~'^. This shows that for 
sources with extreme obscuration, even X-rays are strongly 
attenuated. 



3 X-RAY SPECTRAL ANALYSIS 

The X-ray spectra of galaxies and AGN are very compli- 
cated and as such it is necessary to utilise high signal to 
noise spectra with good spectral resolution and a wide band 
pass in order to correctly determine the source properties. 
These have been provided by XMM-Newton better than ever 
before, so in using data obtained from this observatory, we 
aim to gain new and improved insights into the X-ray spec- 
tral properties of the 12MGS. 



Rush et al. (19961 studied 



the soft X-ray spectra of the 12MGS from ROSAT All-sky 
survey observations, but in doing so will have missed any 
moderate to heavily obscured emission, and furthermore, 
the signal to noise and resolution of ROSAT was incompa- 



rable to that of XMM-Newton. Barcons et al. ( 1995 1 used 
hard X-ray data from HE AO 1 and will therefore have been 
sensitive to obscured sources. However, their data also suf- 
fered from poor signal to noise and low spectral resolution, 
meaning they could only determine the source intensity with 
no meaningful investigation of the absorption or continuum 
slope. 

We develop a systematic approach to the spectral fit- 
ting of our sample in order to determine the intrinsic X-ray 
continuum slope, F, the neutral absorption column density, 
A'h, and the 2-10 keV intrinsic source power, Lx. Charac- 
terising F is necessary for insights into the primary X-ray 
generation process in AGN, thought to be the Comptoni- 



sation of accretion disk photons by a hot corona (Sunyaev 
|fc Titarchu"k||T980} [Haardt fc Maraschi|[T99T| . This intrin- 
sic source continuum parameter is also useful for comparing 
different source populations, such as Seyfert Is and 2s, to de- 
termine if their central engines are the same, as predicted by 
unified schemes. Determining the A^h distribution in AGN 
is a key ingredient to XRB synthesis models, and measur- 
ing the intrinsic power allows for comparison to the power 
at other wavelengths. We will discuss most of these issues 
in a subsequent publication however, and we focus on an 
accurate determination of these parameters here. 
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Figure 4. The flux suppression factors for the spherical model 
(top) and the torus model (bottom) for the observed frame 2- 
10 keV, 10-40 keV, 20-100 keV (INTEGRAL) and 14-195 keV 
(Swift/BAT) bands as a function of Afjj at redshifts of 0, 1, 2 and 
3. The flux suppression factor is defined as the unabsorbed flux 
in that band divided by the observed flux, for F = 2. 



3.1 Sample selection and data reduction 

The sample of observations presented here consists of 
all EPlC-pn XMM-Newton observations of the extended 
12MGS ( [Rush et al^[T993l ) that were public in the XMM- 
Newton data archive as of December 2008, including sources 



which were observed serendipitously. In the case of multi- 
ple observations of the same source, we use the longest ob- 
servation available. The data were reduced homogeneously 
and systematically using the method described by |Nandra| 
et al. (20071 using SAS v7.0 tasks. This involves extraction of 



source events from a circular region, radius 35 arcsecs, cen- 
tred on the source position, and extraction of background 
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events from rectangular regions, ensuring no contamination 
from interloping sources. We search for pileup effects in the 
observations by comparing the distribution of single and 
double pixel events, which are sensitive to these effects, to 
model distributions. Pileup effects are reduced if and when 



detected, by using only CCD pattern (single) events (BaL 



let 



1999). Background flares are filtered out by calculating 
the level of the background count rate at which the ex- 
cess variance is determined to be zero. A background count 
rate of twice this level is used to filter the observation of 
flares. The subsequent background-subtracted spectra ex- 
tracted are then binned with a minimum of 20 counts per 
bin, appropriate for fitting. We include for analysis only 
those sources for which a meaningful spectrum has been pro- 
duced by this process. This excludes those observations for 
which the flaring of the background has rendered the obser- 
vation useless, where the source has not been detected above 
the background, or the detection is at low significance. We 
determine a meaningful spectrum to be one where the nor- 
malisation of a simple power-law model can be constrained 
to be greater than zero (Ax^ > 2.71) with a fit to the 0.2-10 
keV spectrum. After these considerations, the sample con- 
sists of 126 observations. The names, positions and redshifts 
of these galaxies along with the details of the XMM-Newton 
observations are given in Table [2] This set of XMM-Newton 
observations is quite varied, with exposures varying from 
hundreds of seconds to tens of kilo-seconds (Fig. [5|, with 
differing observation modes and filters used. We do however, 
treat each spectral fit the same, regardless of the number and 
distribution of counts in the spectrum (see Fig. [5] for distri- 
bution of counts in these observations and the observed flux 
distribution) . 



3.2 Spectral fitting method: 2.5-10 lieV 

We use XSPEC vll.3 to carry out the spectral fitting of 
the 12 micron selected galaxies. As each spectrum has been 
grouped with a minimum of 20 counts per bin, the errors on 
the data points can be considered as Gaussian, and therefore 
fitting can be used to find the best fit parameters for a 
given model. 

Factors complicating the determination of the intrinsic 
parameters and the A'^h are the soft excess components, ubiq- 
uitous in the soft X-ray spectra of AGN, which may arise 
from the accretion disk, thermal emission, scattered intrinsic 
emission and/or unresolved emission lines. For this reason, 
where possible, we determine the intrinsic source parameters 
from a fit to the 2.5-10 keV spectra for these galaxies, where 
the soft excess components are negligible. The following is 
a description of the spectral fitting methods we use. 

We begin by fitting the 2.5-10 keV spectrum with a 
simple power-law model (power law in xspec). We then sys- 
tematically include further model components in order to fit 
the more complex features that the spectrum presents. We 
note the change in the x'^ value of the fit after the addition 
of each added component and the change in the number of 
degrees of freedom (DOF). For a model combination to be 
allowed, we stipulate that A^^^ 4.0 for each DOF lost in 
adding the additional component (based loosely on 95% con- 
fidence levels quantified using an F-test). The combination 
of model components which meets these criteria and has the 
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Figure 5. X-ray observation statistics for the sample, with expo- 
sure time distribution (top), source count distribution (middle) 
and observed source 2-10 keV flux (bottom) 



smallest reduced x^ (Xr ~ X^/DOF) is presented as our best 
systematic fit for that spectrum. 

We use the following models in combination to acquire 
a fit to the 2.5-10 keV spectrum: 



• Neutral absorption model, zwabs ( Morrison & McCam- 
|mon|1983| l, for absorption intrinsic to the source. 

• Neutral refiection from Compton thick material, pexmon 



( Nandra et al. 2007 1, which is derived from the pexrav model 



of Magdziarz & Zdziarski ( 1995 1 and includes self consistent 



Fe Ka, with its Compton shoulder, and Fe K/3 line emission 
( George fc Fabian||1991| |Matt|2002[ ). 
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• Our new model describing a heavily absorbed power- 
law, trans, based on Monte-Carlo calculations, with a spher- 
ical geometry, which take into account Compton scattering 
and includes line emission (as described in Section [2|. We 
use this model in addition to the zwabs model as it is more 
likely to pick out a heavily obscured source due to the in- 
clusion of line emission. 

• Additional line emission, modelled by a gaussian: 



The pexmon and trans models include Fe Ka line emis- 
sion, produced by cold distant matter, and therefore, if a 
narrow emission feature exists at 6.4 keV, these models 
should be preferentially included in the fit. We also check 
to see if any emission is present from ionised iron by adding 
narrow gaussians at 6.7 and 6.96 which correspond to Ka 
emission from He-like and H-like iron respectively. We also 
add a third gaussian component with energy constrained to 
be 6.2 <Eii„e < 6.6 keV and with a free width parameter, to 
check for broad Fe Ka. This represents AE=0.2 keV from 
the rest-frame energy of neutral Fe Kq so as to avoid fit- 
ting ionised Fe lines at higher energy. In reality, broad Fe 
Ka lines have been foimd with a lower energy than the en- 
ergy we permit (e.g. [Nandra et al.|[2007[ ), but we use this 
conservative range as broad Fe Kq lines are not the focus 
of this investigation. We keep these gaussian components if 
the normalisation of each line is constrained to be greater 
than zero (Ax^ > 2.71). 

We keep the power-law indices for the pexmon and trans 
components equal to the primary power-law for each fit as 
it is presumed that each component in this band has the 
same origin. For spectra where the power-law index, F, is 
badly constrained we fix the power-law index to F = 1.9, 



the canonical value for Seyfert galaxies (Nandra & Pounds 



1994), which make up the majority of our sample, but also 



roughly appropriate for galaxies ( Ptak et al.|1999 1. We de- 
fine fits where F cannot be constrained within AF 5C 0.7 to 
be unconstrained as such an uncertainty introduces a large 
overall uncertainty into the fit and other fit parameters as 
well as the source luminosity. 

We fit the 2.5-10 keV spectra initially as this band is 
least affected by galactic and other complex soft X-ray emis- 
sion, and is therefore best suited for determining the intrinsic 
parameters of the nuclear source. Furthermore, the higher 
signal-to-noise ratio at soft energies often results in spectral 
fits being driven by the soft features. For example, strong 
soft emission may cause the power-law index to steepen and 
as a result the measured A^h may increase to compensate 
for this. However, the lower signal-to-noise ratio in the hard 
band can sometimes lead to unconstrained fits. We do not 
use fits where the normalisation of the power-law is not con- 
strained to be greater than zero and instead use the full band 
fit method, described below. For fits that produce a low F, 
a hidden transmission or reflection component is likely to 
be present and not recognised by the 2.5-10 keV fit. Where 
F < 1.4, this presents a scenario where the photon index 
lies greater than 3-cr away from the mean measured by |Nan-| 
dra & Pounds (19941. In this case we also refer to the full 



3.3 Spectral fitting metiiod: 0.2-10 keV 

Ignoring data below 2.5 keV leaves us insensitive to neutral 
columns less than ~ 10^^ cm^^. For satisfactory hard band 
fits, we then freeze the hard band fit and continue fitting 
the 0.2-10 keV spectrum with these parameters frozen, with 
the exception of the Nu parameter. The following models 
are then systematically added, again requiring that Ax^^ 
4.0 for each degree of freedom lost in adding the additional 
component. Again, the combination of model components 
which meets these criteria and has the smallest reduced 
iXr = xVDOF) is adopted as our best systematic fit for the 
0.2-10 spectrum. 

• Neutral absorption (wabs) associated with the Galaxy, 
frozen to the value given in Table[3]and applied to the whole 
spectrum. This component is assumed to be present regard- 
less of the Ax^ • 

• Absorption applied to the whole spectrum and intrinsic 
to the host galaxy, be it neutral or ionised (see below). If 
no neutral absorption is detected in the hard band, we also 
check to see if a neutral absorption component is required 
for the primary power-law in the soft band (zwabs). 



• A warm absorber (cwalS, based on XSTAR (Kallman 
et al.|2004 l. See Nandra et al. (20071 and references therein 
for description) to model cases where ionised absorption fea- 
tures are present in the spectrum. If evidence for neutral 
absorption is present in the hard band fit, we also check 
whether an improved fit results by adding a warm absorber 
on top of this or indeed, replacing this with a warm absorber. 

• A thermal plasma component (apec, [Smith et al.|2001 1. 
The soft X-ray spectra of obscured AGN have been shown 
to consist partly of emission lines photo-ionised by the AGN 
( Guainazzi fc Bianchi|2007 |. Models such as apec are often 
used to account for this emission when these lines are not 
resolved. When using this model, the abundances are ini- 
tially frozen to the solar values, however, once the best fit 
is obtained, we check if thawing the abundances produces a 
better fit. 

• A second power-law, where the power-law index, r2 is 
constrained such that F2 = Fi . This component is used to fit 
any 'scattered power-law' often seen in Seyfert 2 type spectra 



(Turner et al. 19971. It can only be seen when moderate or 



band fit, where a refiection or heavily obscured transmission 
component is often then detected. 



high absorption is occurring to the primary power-law, and 
therefore we only add it if absorption is detected in the hard 
band fit. 

• We also check to see whether F2 ^ Fi produces a better 
fit. 



3.3.1 Modelling the soft-excess seen in Seyfert 1 like 
spectra 

Modelling the soft-excess seen in quasar and Seyfert 1 like 
spectra is complicated, especially as we do not fully under- 
stand the nature of this feature. The models used to fit 
the soft excess are additive (often a blackbody or Comp- 
tonised component), thus accounting for flux in the soft X- 
ray band. However, if an absorption component is used in 
conjunction, adding the soft excess component can lead to 
an over-estimation of the A^h parameter due to 'balancing' 
of emission and absorption. Thus, we apply strict criteria in 
order to include a soft excess model in the spectral flt. 

This soft excess is thought to originate from close to the 



XMM survey of the 12MGS 9 



central engine, and tiius sliould be an excess above tlie in- 
trinsic, unabsorbed power-law spectrum. To test for the soft 
excess, we extrapolate our fit of the 2.5-10 keV spectrum 
into the soft band, removing any absorption components. 
Only then, if there is emission in excess of this extrapolated 
fit in the soft band, do we include soft-excess components in 
the soft band fit. To detect excess emission in the soft band, 
we use a 'soft ratio', defined Eis the data-to-model ratio at 
0.4 keV, done similarly by [Sobolewska fc Done| ( |200 7'). An 
energy of 0.4 keV is used to avoid atomic emission from oxy- 
gen at ~ 0.7 keV. We choose a soft ratio lower cut off of 1.2 
(20%), which is high enough to exclude uncertainties in the 
spectral calibration (see below). We also require this ratio to 
be rising towards softer energies, as emission from thermal 
plasmas can cause an excess at 0.4 keV, but that peaks at 
higher energy. Figure [6] gives two examples where this soft 
ratio is greater than 20%, but one shows a classical soft ex- 
cess which rises from 1 keV to 0.5 keV (3C273), whereas the 
other shows an excess due to the emission from a thermal 
plasma, which falls from 1 keV to 0.5 keV (Messier 100). If 
the soft ratio exceeds 1.2 and the ratio is rising from 1 to 
0.5 keV, we then proceed to include soft excess model com- 
ponents in the spectral fit. The model components we use 
to fit the soft excess are either a black body (bbody) or a 
Comptonised (compTT) component, whichever provides the 
best fit. 



3.3.2 Bad spectral fits 

We define fits which give Xr > 2.0 as 'bad'. Although this 
definition is arbitrary, it is useful for identifying cases where 
the models used are not a good representation of the spec- 
trum. For the 2.5-10 keV fits, only NGC 1068, NGC 1365 
and M82 have bad fits. These sources have notoriously com- 
plex spectra, which are high signal to noise in this case, 
with several emission and absorption features. We do not 
attempt to fit these additional complexities, however, and 
present only the basic hard band parameters for these fits. 
No errors are quoted because the standard formalism does 
not apply unless the model fits the data ( jLampton et al.| 
1976|. 



When we fit the full band spectrum, we sometimes find 
that keeping the hard band parameters fixed will cause a 
bad full band fit, so we then thaw these parameters for an 
improvement. After this procedure, 17/126 fits are still bad. 
These bad fits usually result from high signal-to-noise spec- 
tra which reveal highly structured soft X-ray features which 
are too complex for our models. Ifowever, the nature of the 
soft X-ray emission is not the focus of this investigation, as 
we seek mainly the intrinsic hard X-ray power of the source, 
along with any neutral absorption present. In these cases, 
we do not present the results from the full band fit, rather 
only the hard band fit. These are denoted by an '•"'> in the 
spectral fit table. Furthermore, we note that the calibra- 
tion uncertainties on XMM-Newton are at the level of ~5% 
( Kirsch et al.|2004 ). If the ratio of the data to the model is 
less than this in the soft band, the bad fit may be partially 
due to calibration uncertainties. In this case we present the 
soft band fit results, despite a bad ■ We present the pa- 
rameters F and the power-law normalisation from the hard 
band fits here as these fits had 'good' Xr values and thus 
the uncertainties could be estimated. These are denoted by 





Figure 6. The ratio of the data to the extrapolated 2.5-10 keV 
power-law for 3C273 (top) and Messier 100 (bottom). 3C273 has 
a classical soft excess at up to more than twice the level of the 
extrapolated power-law. Messier 100 has strong emission from a 
thermal plasma, but this peaks at ~ 1 keV. 



(-D) 



in the spectral fit table. 3C120 is an example of a 



source with a bad spectral fit, but a data/model ratio of 
less than 5%, which we show as an example in Fig. [s] Our 
procedure should not result in grossly inaccurate estimates 
of F, Lx or TVh in such cases. 

For cases where the hard band fit is unconstrained (i.e. 
the power-law normalisation is not constrained to be greater 
than zero) or the hard band fit is bad {Xr > 2.0), then 
we cannot use the method described above, whereby the 
hard band fit is extended into the soft band. In these cases, 
we carry out simultaneous full, 0.2-10 keV band fits with 
all model components. Although this is less desirable, it is 
useful in the case where the hard band fit is unsatisfactory 
and for cases where the hard band measures a flat power-law 
index. These are denoted by an in the spectral fit table. 



3.4 Intrinsic power and reflection fraction 

We calculate the intrinsic power emitted in the 2-10 keV 
band for each galaxy based on the absorption corrected flux 
in the 2-10 keV band from the spectral fits to the primary 
power-law. We define the primary power-law as the power- 
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Figure 7. (top) 2-10 keV intrinsic (absorption corrected) lumi- 
nosity plotted against redshift for all objects, (bottom) The lu- 
minosity correction {Lyi{int) / L-x_{obs)) against the intrinsic X-ray 
luminosity. 



law component with the greatest normalisation, i.e. the in- 
trinsic flux at 1 keV. This may include the trans compo- 
nent, if this is present, which will typically model a heav- 
ily obscured hard component visible only at hard energies. 
Standard cosmological parameters are assumed throughout 
{flm = 0.3, A = 0.7, Ho = 70 km s"^ Mpc"^). In the top 
panel of Fig. [7| the intrinsic, absorption corrected 2-10 keV 
luminosity of each galaxy is plotted against its redshift. This 
figure illustrates that the 12MGS is a low redshift galaxy 
sample, with a wide range of X-ray luminosities from 10^* 
ergs~^ through to 10*^ erg s~^ (though with one very lu- 
minous source, 3C273 at ~ 10'*® ergs~^). The lower panel 
shows the size of the correction made to the observed lumi- 
nosity against the intrinsic luminosity, which can be up to a 
factor of ~50. 

For reflection dominated spectra, we derive the intrinsic 
power of the source from the intrinsic power-law required to 
produce the reflection component. The pexmon model calcu- 
lates the reflection from a semi-inflnte slab, subtending 2n 
steradians on the sky. The estimate of the intrinsic power in 
the reflection dominated sources is uncertain as it depends 
on the precise geometry of the reflector. This was shown by 



Murphy & Yaqoob ( 2009 1 who modelled reflection from a 
toroidal distribution of matter, finding that the reflection 
fraction can be underestimated in this way by a factor of 



~ 6. We also investigate the refiection fraction and intrin- 
sic luminosity derived using the pexmon model and compare 
these to those derived from our own torus model. 



3.5 What can we determine using the torus 
spectral model? 



For spectra in which we find reflection is strong (Rp, 



exmon 



> 



pexmon component to the normalisation of the primary 
power-law), we flt the same spectrum with model spectra 
from the torus simulations with two aims. Firstly, as men- 
tioned above, we aim to investigate the reflection fraction 
and intrinsic source luminosity as determined using the torus 
model with respect to the pexmon slab model, and secondly, 
we investigate if it is possible to determine and constrain 
the torus geometric parameters, 9tor and 6i from the 2-10 
keV spectra. For determining the reflection fraction, we fix 
dtor ~ 60°, leaving 9i free to vary. We then determine if 
the inclination angle can be constrained. Following this, we 
fix 9i = 90°, and allow 6tor to vary, to see if this can be 
constrained. 



3.6 X-ray AGN in our sample 

The main goal of this work is to ascertain the properties 
of the primary emission, as well as the absorption charac- 
teristics of our sample. We do this initially for our whole 
sample, and we then go on to do so for AGN only. We de- 
fine a set of unambiguous X-ray AGN as galaxies with an 
observed 2-10 keV luminosity greater than 10^^ ergs"'^ as 
these sources are unambiguously powered by the accretion 
of material onto a super-massive black hole. No local pure 
star-forming galaxy has ever presented a 2-10 keV luminos- 
ity above this limit. One of the most X-ray luminous star 
forming galaxies is NGC 3256, which presents a 2-10 keV 
luminosity of 2.5 xlO** ergs~^, but shows no evidence for 
AGN activity ( |Moran et al.||1999[ ). Therefore, galaxies with 
an observed 2-10 keV luminosity greater than lO**^ ergs~^ 
are very likely to host an AGN. However, this will not se- 
lect all AGN in this sample, so the remaining sources cannot 
be thought of as non-AGN. In this way 60/126=48% of the 
sources in our sample are unambiguous AGN taking into 
account their X-ray luminosities only. 



4 SPECTRAL FITTING RESULTS 

4.1 Spectral fit components and parameters 

From our spectral fitting method, we have determined a best 
fit model for each source, which includes all the absorption 
and soft excess components, giving us an accurate estimate 
of r, A^H and the intrinsic power accounting for the com- 
plexities in the spectra. For 126 spectra, we find that in 62 
of them, we can constrain F on the primary power-law to 
AF < 0.7 (we fix F = 1.9 otherwise). Furthermore, 96 spec- 
tral fits require an absorption component, 39 of which warm 
absorption has been detected. A thermal plasma component 
is ubiquitous in the fits as well, needed to fit 93 spectra. Fe 
Ka emission is detected in 71 spectra, where 36 of these 
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are broad lines (cr > 0.1 keV). A secondary power-law com- 
ponent is required in 44 of our fits, and reflection is also 
required in 44 of the fits. 

The spectral parameters of all fits are listed in Tables [3] 
to|8] Table [3] presents the basic spectral fit parameters such 
as the fit statistics, the parameters of the primary power-law 
and the thermal plasma component, where present, and the 
observed 2-10 keV flux and intrinsic 2-10 keV luminosity. 
Table |4] gives details of any secondary power-law component 
required in the fit; Table [5] gives details of the sources which 
have a soft excess detected and of the models used to fit this 
component; Table |6] presents details of the reflection com- 
ponent needed in each fit, including data obtained from the 
torus model; Table [T] presents the parameters of the warm 
absorber model if ionised absorption has been detected in 
the spectrum; and Table [8] gives details of the line emission 
detected in the iron K regime. Included on all fit parameters 
in which a satisfactory fit has been acquired are the 90% con- 
fidence errors for one interesting parameter (Ax^=2.706). 

A selection of example spectra are presented in Fig. [8] 
These are the spectra of 3C120, NGC 3079, NGC 4593, NGC 
5170, UGC08850 and NGC7479, chosen to display a range 
in signal-to-noise ratio as well as model components. 3C120 
and NGC 4593 have a soft excess present in their spectra as 
determined by the ratio of the data at 0.4 keV to the ex- 
trapolated 2.5-10 keV power-law, and fitted by the compTT 
model. Reflection from neutral material is included in the 
fit, where in the case of NGC 4593, this component accounts 
for the Fe Ka emission seen. 3C120 requires an additional 
gaussian component to model what seems to be broad Fe 
Ka emission. These fits also include a thermal plasma com- 
ponent, however, it is possible that in these cases this model 
component is not representative of real thermal plasma emis- 
sion as this component is not likely to be detected above the 
emission from the AGN, and is rather being used to fit the 
soft complexities of the spectra. A warm absorber has also 
been detected in the spectrum of NGC 4593. 

The spectrum of NGC 3079 is typical of a reflection 
dominated spectrum {iipexmon ~ 16.3) with an intense Fe 
Ka component (940 eV) and is certainly Compton thick. 
UGC 08850 presents a transmission dominated spectrum, 
which is well fitted by our new model treois measuring A''h 
= 3.8 X 10^'^ cm~^, plus a scattered power-law. In the spec- 
trum of NGC 7479 we have detected a hidden transmission 
component (A^h — 2x 10^* cm~^) beneath a reflection domi- 
nated spectrum, though it is not well constrained. The dom- 
inance of the reflection component and high EW of the Fe 
Ka line (768 eV) imply a high column density in support 
the measured A^^h in the transmission component. 

The spectrum of NGC 5170 has a low signal to noise 
ratio and in this case we have had to fix the power-law index 
to 1.9. The only added model component here is low level 
neutral absorption. 

In Fig. [9j we present the 2.5-10 keV spectrum of NGC 
1068, which is an example of where the 2.5-10 keV fit is 
bad. This is a high signal-to-noise spectrum which reveals 
multiple discrete features, possibly due to absorption and 
emission, or even calibration uncertainties. 



4.2 Properties of the primary emission 

The primary power-law is found to have a mean power-law 
index, < F 1.90lo;o7 for the whole sample, and an intrin- 



sic spread of cr = 0.31 



+0.05 
0.05- 



method of Maccacaro et al. 



We use the maximum likelihood 
( 1988j), which takes into account 



the measurement errors on F, to determine both the mean 
and the intrinsic standard deviation of this distribution. 

We also investigate the F distributions for obscured 
and unobscured sources, which is an important test for 
AGN unification. We do this for the entire sample, and 
also for unambiguous AGN with Lx > 10*^ ergs~^. In 
Fig. [To] we show how the distribution of F differs for ob- 
scured sources in these two cases. We again use the max- 



imum likelihood method of Maccacaro et al. ( 1988) to de- 



termine both the mean and the intrinsic standard deviation 
for each of these distributions. We find that for all unob- 
scured sources, < F >= 1.90l[5;JJ^ and a = 0.33to°ol and for 
all obscured sources, < F >= 1.97+o }2 and a = 0.44+°;|!|g. 
This is statistically consistent with the intrinsic power-law 
indices of these two populations being the same. We also 
find that for unobscured AGN, < F >= 1.84l° i4 and 
and for obscured AGN, < F >= 1.90tj^ i^ 



32+"-" 

"••J^-0.08 



and a = 0.38lJ!;|!,?. This is again statistically consistent with 
the intrinsic power-law indices and their intrinsic spreads 
being the same for both obscured and unobscured AGN. A 
Kolmogorov-Smirnov (K-S) test confirms that the probabil- 
ity that these distributions are drawn from different parent 
samples is low (42%) for AGN. 

Furthermore, we find no de pen dence of F on X-ray lu- 
minosity. We show this in Fig. [ll] in the range of 10"^* to 
10*^ ergs~^, where the scatter of F versus X-ray luminosity 
is plotted, as well as the mean values for each logarithmi- 
cally spaced luminosity bins. The vertical error bars show 
one standard deviation spreads in the values. 



4-2.1 Soft-excess modelling 

We find that a total of 15 spectra in our sample require a 
soft-excess model component in their spectral fit. We have 
determined this from two soft ratio criteria, that the soft 
ratio must be > 1.2 and that is must be rising from 1 keV 
to 0.5 keV. 3C120 has a soft excess which does not fit our 
criteria, as it does not have a data/model ratio greater than 
20% at 0.4 keV. However, this source has a clear soft excess 
at 0.7 keV not well modelled by thermal plasma emission. 
We find that a ratio of 3.2 at 0.4 keV is recovered when an 
additional intrinsic absorption component is used. In Fig. |12[ 
we plot the distribution of the soft ratio, which ranges from 
~1.3-7.8 for our sample, though with only one spectrum 



with a ratio greater than 4. Sobolewska & Done ( 2007 1 note 



that this soft ratio ranges from ~2-3 for their sample of 
quasars. The details of these fits are presented in Table |5] 
All but three of the fits prefer a Comptonised component 
over a black-body component. 



4.3 Obscuration and reflection 

4-3.1 Compton thin sources 

We find that 28/126 (22%) of our spectra show no evidence 
for any absorption at all, whereas 77 (61%) spectra present 
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Figure 8. 0.2-10 keV spectra of 3C120, NGC 3079, NGC 4593, NGC 5170, UGCG8850 and NGC7479, shown with best fit 
model from Table [s] 
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Figure 9. 2.5-10 IceV spectrum of NGC 1068 (left), shown with best fit model (right) from Tablejs] This is an example of a spectrum 
where the fit is bad due to spectral features too complex for the models we use. 



Compton thin (neutral) absorption, showing that absorption 
is commonplace in the X-ray spectra of all galaxies. 

4-3.2 Compton thick sources 

As mentioned previously, for 0.2-10 keV X-ray spectra, a 
direct measurement of a column density of Compton thick 



source is difficult due to the extreme suppression of the X- 
ray flux in this band. We have, however, successfully used 
our new model, trans, to directly measure the A^h to be 
greater than 1.5 x 10^* cm~'^ in four sources (NGC 1320, 
NGC 1667, NGC 4968 and NGC 7479). The inclusion of 
a self-consistently modelled iron Ka line is likely to have 
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Figure 10. Top: V distribution of all sources (empty histogram), 
and obscured sources (A'^h > lO'^^ cm~-^, black histogram). Bot- 
tom: r distribution of all unambiguous AGN, defined by Lx > 
10*^ ergs s~^ (empty histogram), and obscured AGN (A^h > lO^'^ 
cm~^, black histogram). The F presented is that from the primary 
power-law only. 



provided greater constraint on the A^h than models which 
include this component separately. 

In our analysis, we call sources with a reflection dom- 
inated spectrum Compton thick, of which there are 12. 
We define a refiection dominated spectrum to be where 



> 



The R„ 



> 8 criterion is arbitrary and 



chosen to be low enough to include known Compton thick 
AGN such as NGC 1068. Figure [13] shows the EW of the 
Fe Ka line against A^h for our sample for reflection domi- 
nated sources (green), mildly refiection dominated sources 
(blue, 1 <R< 8) and sources with weak or no refiection (red, 
R< 1). It can be seen from this that reflection dominated 
sources also have high iron line EWs. We plot only narrow 
line EWs here, which are most likely to originate in cold 
distant obscuring matter, unless no narrow line has been 
measured, in which case we plot the broad line EW, if one 
has been measured 

There are two sources, however, that have high iron 
line EWs which are constrained to be > 300 keV, but have 



Figure 11. Plot of T vs. 2-10 keV intrinsic luminosity. Red crosses 
give the mean T for each luminosity bin. The vertical error bars 
show one standard deviation spreads in the values. 



•z. 




Figure 12. Distribution of the soft ratio, the data-to-model ra- 
tio at 0.4 keV to the extrapolated 2.5-10 keV power-law, for 15 
sources found to have a soft excess. 



a measured Nyi< 10^^ cm~^. Our Monte-Carlo simulations 
have shown that this is not possible for either a spherical or 
toroidal matter distribution, and so we investigate them fur- 
ther. These sources are NGC 3690 and IRASF13349-I-2438. 
NGC 3690 is the merging system ARP 299, which was shown 
by Beppo-SAX to host a Compton thick AGN ( [Delia Ceca| 



et al.,2002), confirming our assertion. We therefore add it 



to our sample of Compton thick AGN. |Longinotti et ah] 
( 2003 1 examine the Fe Ka features in the X-ray spectrum of 
IRASF13349-I-2438, a narrow-line Seyfert 1, finding them to 
be broad and complex. This line is therefore unlikely to orig- 
inate from distant Compton thick obscuring matter. This 
shows that Fe Ka EW diagnostics for obscuration are only 
suitable when we can be sure that the line originates in cold, 
distant matter such as the torus. 

For sources where their Compton thick nature has been 
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Figure 13. The equivalent width of iron Ka line versus directly 
measured Afjj for our sample. The lines drawn show the theo- 
retical relationship between the EW of the Fe Ka line and the 
Affj, calculated from the Monte-Carlo simulations, described in 
section [2] for a power-law with r=2 and solar abundances. The 
solid blue line shows the calcSlations for the spherical distribu- 
tion, the dashed green line shows the calculations from the torus 
distribution for the equatorial line of sight, and a torus opening 
angle of 0tor=6O°. The dotted yellow lines show maximum and 
minimum EWs observed from unobscured sight lines in the torus 
distribution, taken from Figls] where the line of sight A'^h through 
the torus is greater than 10 cm~^. Thus, all truly unobscured 
sources must have an EW less than the maximum shown. Sources 
which are reflection dominated (R>8) are shown in green. Sources 
with 1 <R< 8 are shown in blue for comparison. All grey lines are 
error bars. These are only shown where the EW is consistent with 
being non-zero (plotted as filled diamonds, otherwise plotted as 
open diamonds). Arrows indicate upper limits on the A'^h- 



inferred by the EW and the reflection fraction, we can only 
put a lower limit on the A^h of 1.5 x 10^* cm~'^, but in 
reality, this could be higher and requires data above 10 keV 
to constrain. We find a total of 16 Compton thick sources in 
our sample, the details of which are given in Table [lO] The 
Compton thick fraction for all 126 galaxies is 16/126 (=13 
± 3%), where the errors have been calculated for a binomial 
distribution [op — •\/p(l — p)/N , where p is the proportion 
and N the total number of sources). The Compton thick 
fraction for unambiguous AGN (Lx > lO*^ ergs~^) is 11/60 



(=18 ± 5%). 



4-3.3 Nh distribution and obscured fraction 

In Fig. |14| we present the A^h distribution for the whole 
sample, and over-plotted is the A''h distribution of those 
with Lx > 10*^ ergs"\ The main difference in the two 
distributions is the strong peak of unobscured sources at 
low X-ray luminosities and the greater proportion of heav- 
ily absorbed X-ray luminous AGN. The low luminosity, un- 
absorbed sources may consist of unobscured star forming 
galaxies, LINERs or lower luminosity AGN, but we will in- 
vestigate this further in paper II. The A^h distribution of 
unambiguous AGN shows a smaller unobscured peak, but 
with a second obscured peak at 7Vh = 10^'^~^''. 



Figure 14. distribution of all sources (empty histogram), 
and unambiguous AGN, defined by Lx > 10*'^ ergs s~^ (black 
histogram). The Nn presented is that measured from the primary 
power-law only. 



Fig. [15] plots A^H vs. Lx and also the obscured fraction, 
defined as the ratio of the number of sources with A'h > 10^^ 
cm~^ to the total number of sources in each luminosity bin. 
There is significant variation in the obscured fraction with 
X-ray luminosity, particularly above and below 10'*^ ergs~^ 
where the ratio declines steeply from a peak. 



4.3.4 Reflection dominated sources 

For sources where the X-ray spectrum includes a strong re- 
flection component, we have also fltted our new torus model, 
flrstly to investigate the reflection fraction and intrinsic 
source luminosities as determined using this model with re- 
spect to the pexmon slab model. Secondly, we investigate if it 
is possible to determine the torus geometric parameters, 9tor 
and 9i from the 2-10 keV spectra. In Fig.[T6]we plot the dis- 
tributions of the ratio of the reflection fraction as determined 
by the torus model to the reflection fraction determined by 
the pexmon model and the ratio of the intrinsic luminosities 
as determined by each model. We flnd that for most sources, 
both the reflection fraction and intrinsic source luminosity 
as determined using the torus model is greater than that 
determined by the pexmon model. The means and medians 
of these ratios are 3.08 and 2.70 for J?torus/i?pcxmon and 2.62 
and 1.99 for LXtorus/LXpcxmon. We also flnd that it is not 
possible to constrain the inclination angle or opening angle 
of the torus with these data below 10 keV, as neither pa- 
rameter can be constrained, even with the other parameter 
flxed. 



5 DISCUSSION 

5.1 X-ray transmission in heavily obscured 
sources 

We have made predictions of the equivalent width (EW) 
of the iron line for different model parameters and for two 
matter geometries, a sphere and a torus. The predictions 
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Figure 15. A^h and obscured fraction versus 2-10 keV luminosity 
for all sources. Obscured fraction is calculated in equal logarith- 
mic luminosity bins, as the number of sources with Afjj > 10^^ 
cm"'^ divided by the total number of sources in that bin, where 
more than one source exists per bin. 
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Figure 16. Histograms showing the distribution of (top) the re- 
flection fraction and (bottom) intrinsic 2-10 keV luminosities as 
determined from the torus model with respect to the pexmon 
model. 



we make for the torus are useful for comparing observed 
EWs and making inferences from these about the matter 
geometry and parameters in the observed system. For ex- 
ample, we find that a maximum EW of 150 eV is achievable 
when considering a torus geometry for unobscured sight- 
lines with a typical V (2.0) and solar abundances. To achieve 
greater EWs than this, the line of sight must be obscured 
with NyO 10^^ cm"^. Therefore, if an EW higher than 
this is measured for an unabsorbed spectrum, it may cast 
doubts on the reliability of the measured line of sight A^h. 
We showed this for NGC 3690, which shows an Fe Ka EW of 



913_34g eV, but a measured A^h of only 6.14_g 



xlO^ 



which is clearly at odds with our model predictions. Data 
above 10 keV from Beppo-SAX shows, however, that this 
source is indeed heavily obscured with A^h> 1.5 x lO^^cm"'^ 
( [Delia Ceca et aT]|2002[ ), in full support of our model pre- 
dictions. Our Fe Ka EW predictions also agree well with 



predictions of the same nature as made by 'Ghiscllin i et ah] 
( |1994| ), [Ikeda et al.| ( [20091 ) and [MurphylTyaqoob, (2009| . 

We also present flux suppression factors for both spher- 
ical and toroidal geometries as a function of A^h for differ- 
ent bands and redshifts. These are particularly relevant for 
AGN surveys that use the X-ray band as they indicate the 
effectiveness of each band at detecting heavily obscured sys- 
tems. It is known that even the hardest X-ray energies are 
affected by extreme obscuration. Here we quantify the flux 
suppression in these bands, finding that locally, flux suppres- 
sion reaches a factor of 10 at A^h=10^^ cm~^ for the > 10 
keV bands. This may explain the fewer than expected local 
Compton thick AGN detected in hard X-ray surveys con- 
ducted by Swift/BAT and INTEGRAL. For the > 10 keV 
X-ray bands, the flux suppression factor is 10 at ~ 10^^ cm"'^ 
for all redshifts. This is because at these energies, Compton 
scattering is the dominant interaction and source of attenu- 
ation in the spectrum. The cross section of this interaction 
is roughly constant as a function of energy and therefore 
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the effect of redshift is negligible. Flux suppression in the 
2-10 keV band reaches 10 at 6 x 10^^ cm~^ but this im- 
proves as redshift increases, becoming ~ 10^^ cm~^ at z=3. 
For the torus geometry, the flux suppression factors behave 
similarly to the spherical geometry up to ~ 4 x 10^^ cm~^. 
Here, the flux suppression factors level off due to reflected 
flux from the inner part of the torus, which remains constant 
regardless of the line of sight A^h- In all cases though, an ab- 
sorbing column of A'^h^IO^^ leaves all X-ray bands severely 
insensitive to detection of X-ray sources. 

Using model spectral data from the torus simulations, 
we also investigated the difference in the reflection fractions 
and intrinsic luminosities as determined by using these new 
torus data, to that derived from fits with the slab model, 
pexmon. Murphy & Yaqoob ( 2009 1 find that the reflection 
fraction is underestimated by a factor of ~ 6 when they 
compare their torus data to slab data. We also find that the 
reflection fraction is systematically higher when using the 
torus model data, by a factor of 3 on average. Though this 
is not as high as the underestimation presented by |Mur-| 
phy & Yaqoob (20091, the discrepancy is probably due to 



19951. If a systematic difference is found in these distribu- 



differing torus geometries, as they use a cylindrical torus 
geometry which gives varying A^h as a function of inclina- 
tion angle, whereas we use a spherical torus with constant 
A'h. We also flnd though, that the intrinsic luminosity de- 
rived from the slab models are underestimated with respect 
to the torus model by a factor of 2-3. Additionally, using 
the torus model, we find that we are unable to constrain the 
torus geometrical parameters, the opening angle or the in- 
clination angle using the XMM-Newton spectral data. |Ikeda| 
( 2009 1 also find that it wasn't possible to constrain 



et al. 



these parameters with their model, even with spectral data 
above 10 keV with Suzaku 



5.2 Properties of the primary X-ray emission 
from 12 micron selected galaxies 

We have presented the properties of the primary X-ray emis- 
sion of our sample of 12 micron selected galaxies. As our 
sample is a heterogeneous mix of galaxies, these results re- 
flect the properties of all X-ray sources. We will present re- 
sults on individual optical classes in paper II. We do however 
deflne a subset of unambiguous AGN using their observed 
2-10 keV luminosity though this will of course not include 
AGN with low X-ray luminosity which may be due to sig- 
nificant absorption or being of intrinsic low luminosity. 
We find that the mean power-law index is < F 



1.90 



+0.05 
-0.07 

0.31 



for the whole sample, with an intrinsic spread of 
The mean is consistent with previous works 



+0.05 



on predominantly X-ray selected samples. For example |Nan-| 
dra & Founds! (119941) who find < F 1.95 ± 0.05 and 



, = 0.15 ± 0.04 for their sample of Seyfert galaxies ob- 
served with Ginga. We do find a larger intrinsic spread in 



the indices though than Nandra & Pounds (19941. Our re- 



sult is also consistent with the results from harder X-ray se- 
lected samp les (eg. 1.93 for afl Sey ferts in the INTEGRAL 
catalogue of Beckmann et al.|2009 l. 

Comparing the distribution of power-law indices for ob- 
scured sources and unobscured sources is a key test of the 
unification scheme, which states that the central engine is 
the same for both, the difference being the orientation of the 



tions, it would suggest that this is not simply the case. Com- 
paring the average X-ray spectra of Seyfert galaxies with 
Ginga and OSSE data, [Zdziarski et al., ( ,1995j originally 
found a systematic difference between the F distributions 
of Seyfert Is and Seyfert 2s, however, more recent works 
have tended to find that the distributions are indeed con- 
sistent with each other, generally due to the inclusion of 
more complex spectral models, such as Compton reflection 



components (e.g. Dadina 20081. We investigated these dis- 
tributions for obscured and unobscured sources, as well as 
obscured and unobscured unambiguous AGN. Our analysis 
showed that the means of these distributions in both cases 
are not significantly different (the means lie within one sigma 
of each other) and a K-S test supports the same conclusion, 
implying that the X-ray generation mechanism in both ob- 
scured and unobscured AGN is the same, lending support to 
AGN unification. The fact that the intrinsic dispersions are 



similar also supports this conclusion. Beckmann et al. ( 2009 1 



carried out the same analysis for obscured and unobscured 
sources on INTEGRAL data, which should be less affected 
by absorption effects, and also find no significant differences 
between their obscured and unobscured sources. Despite all 
of this, however, we do notice that the F distribution of ob- 
scured AGN is skewed towards low F values. We discount 
the possibility of unidentified refiection in these sources, but 
results presented by Malizia et al. (2003) suggest an alterna- 
tive cause. They also find a harder distribution of F values 
for Seyfert 2s with respect to Seyfert Is from their analysis of 
Beppo-SAX data, but reconcile this difference by including 
more complex absorption models, such as partial covering 
scenarios, which also improve their fits. There have been 
several instances where absorption in AGN has been shown 
to be more complex than a single smooth obscurer, for ex- 
ample NGC 1365, which shows rapid variability in the line 
of sight obscuratio n |Risaliti et al.[2009] ) and further exam- 
ples by Ricci et al. (20101. As we do not use partial covering 



models, this may be the reason for the skewed distribution 
towards low F values in our obscured AGN. 



We have also investigated any correlation between the 



power-law index and the intrinsic X-ray luminosity (Fig. 11 1 



Our results don't support a correlation between these two 
quantities however. This was also found by [George et al.| 
(2000) in their sample of local AGN using ASGA data. At 
higher redshifts however, there have been reports of a corre- 
lation between F and Lx . For example Dai et al. ( 2004 1 find 



an anti-correlation between the two for their mini-survey of 
1.7 < z < 4 sources observed with XMM-Newton and Chan- 



dra in the luminosity range of 10 



ergs . This finding 



was later supported by Saez et al. ( 2008 1 for the Chandra 
deep fields who use a larger range in luminosity (10*^"*^ 
ergs"^) and redshift (0.1-4) with 173 AGN. Our luminosity 
range is io^*~^* ergs~^, which is somewhat lower than the 
higher redsh ift studies, though the local sample of [George] 



system to the observer (Antonucci 1993 Urry & Padovani 



et al. (20001 ranges from 10 ~ ergs" , with no correla- 
tion. This may suggest an evolutionary effect, or a selection 
effect in the X-ray samples. Our sample does not suffer from 
such effects, as it is MIR selected. 
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5.3 Properties of X-ray obscuration in 12 micron 
selected galaxies 

5.3.1 Compton thick sources 

We find that 16/126 (=13 ± 3%) of the whole sample of 
sources studied in this sample are Compton thick. For X- 
ray bright AGN (Lx > lO"^ ergs"^), there are 11/60 (=18 
± 5%) Compton thick sources. This fraction is important 
for synthesis models of the XRB and the accretion history 
of the universe. For X-ray samples selected above 10 keV, 
the Compton thick fraction has been constrained to be less 



than that predicted by the XRB. Beckmann et al. (20091 



found that only 4 ± 2% of the AGN in their sample are 
Compton thick, whereas the XRB predicts a fraction close 
to 15% at the corresponding limiting flux. |Treister et al.| 
( 2009 1 claim that the observed differences between observa- 



tion and models is due to degeneracies in the parameters of 
the models, and that only direct observations of Compton 
thick AGN can determine their relative number. Howeverm 
when accounting for the biases present in hard X-ray se- 



lection, Burlon et al. (20101 recover a 20% Compton thick 



fraction in the Swift/BAT survey, which is consistent with 
the XRB and our results. 

The Compton thick fraction that we measure in 12 mi- 
cron selected sources is almost S-a greater than that mea- 
sured in hard X-ray selected sources, suggesting that mid- 
IR selection is more efficient than hard X-ray selection in 
finding Compton thick AGN. This is explained by the sup- 
pression of X-ray flux in the >10 keV bands at high column 
densities, which we have shown to be a factor >10 at A^h 



= 10^' 
suits. 



cm using our Monte-Carlo X-ray transmission re- 



The XRB synthesis models of |Gilh et aT] ( |2007| ) allow 
the prediction of the number counts expected from an X- 
ray survey given the set of parameters of that survejQ As 
this sample is not X-ray flux limited, the comparison is not 
straightforward. However, the Compton thick fraction ex- 
pected for a survey with our redshift range (0 < 2 < 0.1) 
and 2-10 keV luminosity range (10*^ < Lx < lO"", a lower 
limit of Lx > as we are only considering AGN) lies 

between ~2% and ~20% for X-ray limiting fluxes between 
~ 10"^^ and ~ 10"" ergs cm ^ s ^, in between which our 
'effective' X-ray flux limit lies (see Fig. [sj. Thus our results 
would seem roughly consistent with the XRB models. 

From the four Compton thick AGN which we manage 
to measure the A'h directly with trans, in three of these 
cases, a prominent iron line aids in the determination of the 
A'h. However in the fourth source, NGC 1667, the spectrum 
does not present an prominent iron line. Here the model is 
required to flt a hard excess seen in spectrum, but the Ah 
is not well constrained (Ah = 27lt25i'')- This is a low sig- 
nal to noise spectrum, and should be considered a marginal 
Compton thick source. This technique of fltting hard ex- 
Brightman fc Nandra| ([2008 ) , where we show that 



it can be an effective method of identifying heavily obscured 
nuclei (e.g. the Chandra observation of NGC 4501). In the 
absence of high signal to noise data around the Fe Ka line, 
or data above 10 keV, this method is the next best thing 



An 



online 



calculator 



be 



found 



liere; 



for picking out Compton thick candidates from X-ray spec- 
tral data alone. This may then motivate longer exposures to 
gain higher signal to noise data, or observations of broader 
banded telescopes such as Suzaku. 

We also infer the Compton thick nature of 12 sources 
from their reflection dominated spectra, which we deflne as 
being those with Rpexmon > 8. Though this definition is 
arbitrary, it is well supported in most sources by high Fe 
Ka EWs or a directly measured Ah> 10^* cm~^. Lowering 
this criteron to, e.g. Rpexmon > 1 would include two further 
sources, NGC 4631 and NGC 4666, without prior indications 
of being heavily absorbed. However, these sources don't have 
well constrained Fe Ka EWs and hence don't support a 
lower Rpexmon Criterion. If a higher value was chosen, e.g. 



Rv 



> 20, this would exclude two sources, NGC 1068 



and 2MASXJ15504152. However, NGC 1068 is weU known 
to be Compton thick, and as such motivated the criterion in 
the first place. We therefore conclude that the Rpexmon > 8 
criterion is best suited to selecting Compton thick AGN. 

5.3.2 Nu distribuUon and obscured fraction 

The Ah distribution of AGN has a peak at low column 
densities, which consists of all unobscured sources. For ob- 
scured sources, the distribution peaks in the Ah = lO^^"^** 
cm"^ bin. We find that 29/60 (=48 ± 6%) of X-ray lumi- 
nous AGN are heavily obscured (Ah ^ lO^'^ cm~^). Recent 
studies into the distribution of absorbing columns, such as 



that by Tueller et al. (20081, reveal that for their sample of 
Swtft/BAT 14-195 keV selected AGN, 31 ± 5% are heavily 
obscured, which is a greater than 3-cr difference. Their sam- 
ple does however extend to higher X-ray luminosities, which 
will include a greater fraction of unobscured sources due to 
the dependence they find of obscuration on luminosity. We 
also intend to investigate the Ah distribution for the AGN 
in our sample, but using optical line emission to define the 
AGN in paper II. It is likely here that AGN definition based 
on X-ray luminosity alone will bias a sample. Those sources 
which do not meet the luminosity criterion, but have Ah 
> 10^^ cm"^ are also likely to be AGN. 

Our analysis of the dependence of obscuration on lumi- 
nosity shows that obscuration depends heavily on the intrin- 
sic power of the X-ray source. The data show that obscura- 
tion peaks at Lx = 10'*^ ergs~^, and declines steeply both 
above and below this. Several previous studies have reported 
the decline in the obscured fraction above lO*'^ ergs~^ in 
X-ray selected samples (Ued a et al.| 2003'; 'La Franca et al 



2005 



lAkylaset aLp0 06 Tuo Uer et aL|20 08 Bec kmann et al 



which we confirm here in our MIR selected sample. 
These results are in support of the 'receding torus' model 



of Lawrence (19911 and Simpson (20051 which describe the 



decrease in the covering fraction of the torus in AGN with 
source luminosity, due to sublimation of dust. The decline 
that we show below 10"*^ ergs~^ may be caused in part 
by the inclusion of non-AGN powered sources which emit 
at these lower luminosities. However, recently, [Burlon et al.| 
( [2010 ') find compelling evidence for a decrease in the ob- 
scured fraction at low X-ray luminosities from analysis of 
the X-ray luminosity functions of obscured and unobscured 



http://www.bo.astro.lt/~gilIi/counts.htmI 



AGN from the Swift/BAT survey. Akylas & Georgantopou- 
|los| ( |2009[ ) also evidence for a lower incidence of obscuration 
at low X-ray luminosities in their optically selected sample. 
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Both sets of authors suggest this may be due to the suste- 
nance of torus by accretion driven winds modelled by lElitzur] 
& Shlosman ( 2006!). Our result must be confirmed though 
with a non-X-ray luminosity dependent method of identify- 
ing AGNs, which we intend to do in paper If using optical 
AGN classifications. 



6 CONCLUSIONS 

In summary, we have presented an X-ray spectral analy- 
sis of the 126 galaxies with XMM-Newton coverage in the 
12MGS, which is a relatively unbiased and representative 
selection of galaxies. We have done so with the help of new 
X-ray radiative transfer calculations, based on Monte-Carlo 
simulations which take into account photo-electric absorp- 
tion, Compton scattering and includes fluorescent iron line 
emission. We have compiled new X-ray spectral models from 
these data, trans and torus which we have presented here. 
These table models are available publicly on the web at 
[http : / / astro . ic. ac. u k/mbrightmcm/home 
Our main results from these models have been: 

• unobscured AGN with a presumed toroidal distribution 
of matter around it can achieve a maximum iron Ka equiv- 
alent width of ~ 150 eV, a useful limit to consider when 
assessing if a source is truly unobscured or not. From our 
calculations, we show that in order to exceed this EW, the 
line of sight must be obscured with A^h> 10^^ cm~^ (e.g. 
NGC 3690). 

• for AfH=10^^ cm"^, the flux seen in the 10-40 keV, 20- 
100 keV {INTEGRAL), and 14-195 keV (Swift/BAT) bands 
is only 10% of the intrinsic flux, which is important for con- 
sidering the biases present against hard X-ray selected, heav- 
ily obscured AGN. At this A^h, 10% or less of the intrinsic 
flux is seen in any X-ray band at all redshifts, which reveals 
the bias present in all X-ray surveys against such heavily 
obscured systems. 

• using spectral models based on slab geometries such as 
pexrav or pexmon will underestimate the reflection fraction 
and intrinsic luminosities with respect to toroidal geometries 
by a factor of 2-3, leading to underestimation of the intrinsic 
luminosity of the source, as also found by |Murphy fc Yaqoob| 
(20091). 



By combining the positive attributes of 12 micron selec- 
tion with the penetrative power of X-rays, we have been able 
to determine the intrinsic X-ray properties of galaxies such 
as source power, power-law index, and in addition we have 
assessed the absorption occurring in the sources. We have 
improved upon previous works on this sample with better 
signal to noise, spectral resolution and larger band pass by 
using XMM-Newton data. We have developed a detailed sys- 
tematic approach to the fitting of the spectral data, with the 
aim of determining the intrinsic power, F and A^h as accu- 
rately as possible. We have paid particular attention to the 
uncovering of Compton thick sources using our new model, 
trams and the pexmon model. The main conclusions from 
this study have been: 

• we find the mean primary power-law index to be < 



F 1.90 



+0.05 



for our sample with an intrinsic spread of 



a — 0.31_Q Qg. The index is consistent with previous works. 



though we find a larger intrinsic dispersion in our results 
than previously reported. 

• we also find that the mean power-law index for obscured 
and unobscured sources is consistent with the two source 
populations being the same, supporting unified schemes. 

• we find no dependence of F on X-ray luminosity as has 
been previously reported for local AGN. 

• we find an 18% Compton thick fraction for the X-ray 
AGN defined here, which is higher than the hard X-ray se- 
lected samples, but roughly consistent with XRB predic- 
tions. 

• the obscured fraction for our sample is a strong function 
of X-ray luminosity, peaking at a luminosity of ~ ]^q42-43 
ergs~^. The decline in the obscured fraction at high lumi- 
nosities, where these sources are unambiguously AGN, is a 
confirmation of previous works. The decline at lower lumi- 
nosities has also been suggested but needs further explana- 
tion accounting for optical spectroscopic classifications. 
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Table 4. Parameters of the second power-law component, if it is present. Column (1) Galaxy name; Column (2) the neutral column 
density measured in the secondary power-law (10^^ cm~^); Column (3) the power-law index of the secondary power-law; Column (4) 
the ratio of the normalisation of the second power-law to the primary power-law. 
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Table 5. Parameters of the soft excess component, if it is present. Column (1) Galaxy name; Column (2) the ratio of the data at 0.4 
keV to the extrapolated 2.5-10 keV power-law. * 3C120 has a soft excess which does not fit our criteria, as it does not have a data/model 
ratio greater than 20% at 0.4 keV. However, this source has a clear soft excess at 0.7 kcV not well modelled by thermal plasma emission. 
We find that a ratio of 3.2 at 0.4 keV is recovered when an additional intrinsic absorption component is used; Column (3) the model 

used; Column (4) the input soft photon (Wien) temperature (compTT model only, keV); Column (5) the temperature (keV); Column (6) 
the plasma optical depth (compTT model only); Column (7) the ratio of the normalisation of this component to the normalisation of the 
primary power-law. 
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0.487 tV^l 


NGC4253 


1.62 


compTT 


0.08 


34.1 




0.1 


0.17 


3C273 


2.1 


compTT 


0.06 


66.2 




0.03 


0.14 


NGC4593 


1.64 


compTT 


00743 +0 0257 
U.UUf^O _o. 00486 


2.06 


+0.85 
-0.06 


1- +0.27 
•^■1' -0.56 


590 ttll 


NGC4725 


3.83 


bbody 




0.11 


+0.02 
-0.02 




0676 +0 0679 
u.uom -0.0087 


NGC5548 


1.61 


compTT 


0.02 


26.5 




0.27 


4.44 


UGC09944 


7.8 


bbody 




0.10 


+0.01 
-0.01 




136 

u.iao _(, 014 


MRK0509 


2.87 


compTT 


0.07 


44.6 




0.11 


0.25 


NGC7213 


1.3 


compTT 


146 +°°" 
U.ITO _o.009 


44.9 


+ 12.1 
-42.9 


0.0100 torn 


0141 

-0.0026 


NGC7469 


2.15 


compTT 


0.02 


40.8 




0.07 


8.78 



30 M. Brightman and K. Nandra 



Table 6. Parameters of the reflection component, if it is present. Column (1) Galaxy name; Column (2) power-law index of the reflection 
component; Column (3) the ratio of the normalisation of the pexmon component to the normalisation of the primary power-law; Column (4) 
the ratio of the normalisation of the torus component to the normalisation of the primary power-law. Presented only where Rpexmon > 1; 
Column (5) The intrinsic 2-10 keV luminosity as determined from the pexmon component; Column (6) The intrinsic 2-10 keV luminosity 
as determined from the torus component. 



Name 


Pa 


^pexmon 




LffX (pexmon) 


LffX (torus) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


IVlrxiA-UooO 


O HQ 

z.Uo 


n Qnc +0.115 
0.308 -0.116 








NGC0424 


= ri 


Q 79 +218. 
-9.02 


24.9 


42.5 


43.1 




— i 1 


n IQq +0.105 








NGC0513 


= Pi 


n r»r-r- +1.09 

0.955 +1^^ 




" 




JNGClOoo 


= Pi 


10.3 


31.0 


42.2 


42.4 


ARP1 1 R 


— Pi 
^1 


501 +0-189 
U.OUl -0.213 








NnCI 1 94 


— r. 

— A 1 


4q Q +45.8 
^y-o -46.9 


167.0 


42.3 


43.0 




= Pi 


An A +44.9 
40-4 I39.0 


109.0 


42.7 


42.3 




l.gt 


145 +1'^° 
-■■^O -147 


615.0 


40.9 


41.7 




— i 1 


n it;fl +0.116 








MKivUDio 


= i 1 


n -iTQ +0.507 
u.o/» -0.381 










— 1 


47 +5.24 
-6.11 


11.4 


42.4 


42.4 




— Fi 

— A 1 


441 +0-209 
U-441 -0.168 








iVlx\iVU ( U'i 


— i 1 


n 98"^ +0.708 
-0.290 








NOP9QQ9 


— Fi 

— i 1 


n 220 +°"*'' 
u.zzu -0.109 








iwrpcoTpDriai 

iVlH/OOJ-i-Ji- VUO -L 


— Fi 

— 1 


n '^09 +0-105 
-0.120 










l.gt 


^(^ +7.8 
10-5 -9.7 


30.9 


40.9 


41.2 




— i 1 


n 1 Qn +1.22 
0.139 _m4i 








■\T O O O O "7 


= Ti 


t- 01 +3.02 


12.4 


41.6 


42.0 


i> Kjy^'-i.KJ'.j 1. 


— Fi 

— A 1 


, c:o +1.04 
1.00 _i 55 


9.0 


40.9 


41.5 


Nr;p4i ^1 


— Fi 

— J- 1 


U.(C50 _o.306 










9 OS 
z.uo 


„ -10 +0.082 
0-413 _o.077 








1 Vi J_J kJ kJ A J_J X l/L/ O O 


l.gt 


17 6 +24-2 

I'-" -17.6 


17.5 


41.7 


41.7 


iVlil/ooiiljrvUOo 


— I 1 


n o7r, +0.703 
0-»'^ -0.571 








IN V-TVy^oyo 


— Pi 

— 1 


888 +"-2-'" 










— A 1 


1 16 +2-01 

-1.16 


0.8 


38.9 


38.8 




— ^1 


2 02 +1-58 


8.9 


40.5 


41.1 


NGC5194 


= Pi 


ri +31.6 

51.3 Zli 1 


177.0 


40.4 


40.9 


ESO383-G035 


= Pi 


0.212 +0-06| 

— U.UDo 








1C4329A 


= Pi 


274 +° °^® 
"-z't -0.069 








NGC5506 


= Pi 


270 +°-^'"' 
U.Z(U -0.153 








NGC5548 


1.64 


162 +0-134 
U.IOZ -0.162 








2MASXJ15504152 


1.9t 


10 7 +69.6 
l'^-' -8.9 


56.7 


43.0 


43.6 


NGC6552 


= ri 


1 1 2 +*'^ 
11^ -48 


1020.0 


42.9 


43.6 


NGC6890 


1.9t 


= 9 1 +43.2 
°^-l -47.0 


76.6 


42.2 


42.3 


MRK0509 


= Pi 


309 +° '"^2 
u.oua -0.071 








NGC7172 


= Pi 


255 +°-l°2 
U.zoo -0.256 








3C445 


= Pi 


317 +2-'^'l 
^••^1' -0.152 








NGC7314 


= Pi 


0.131 








NGC7469 


= Pi 


0.698 1°;°^^ 








NGC7479 


= Pi 


293 +3.71 
U.zao -0.295 








ESO148-IG002 


= ri 


41 6 +32.6 
tl-D -43.1 


45.3 


43.2 


43.2 


NGC7582 


1.49 


107 


400.0 


42.6 


42.5 


NGC7674 


1 91 

^•^^ -0.37 


143 +120 
1*'^ -148 


185.0 


43.6 


43.7 



XMM survey of the 12MGS 31 



Table 7. X-ray spectral fitting data. Parameters of the warm absorber, if it is present. Column (1) Galaxy name; Column (2) the 
logarithm of the ionised column density (cm~^); Column (3) the logarithm of the ionisation parameter, ^. 



Name logioA^H logioC 

_a) (2) (3) 

3.88 

+0.8 q OR +0.74 
-0.2 -0.25 

+0.4 r, +0.19 
-0.6 -0.21 
+0-3 A QQ +1.28 
-0.2 ^-"^ -0.00 

+°-3 .0 fiSf) 0.40 
-0.1 U.DOU _{) 55 

+ 0.3 r, r,4f, +0.738 

-0.4 u-^'i" -0.00 

+ 0.2 r, 71 +0.21 

-0.1 '^■'^ -0.23 

+0.2 PI ooA +0.226 

-0.4 U-^O^ -0.00 
+0.0 n on 0.25 

-0.1 -J-yu -0.10 

1°:^ -0.161 

+0.3 o nn +0-32 
-0.7 -0.19 

+0.0 n 1 in +0.082 

-0.3 u-J^J^o -0.028 
+0.0 n fi4 +0.04 

-0.1 -•J-O'* -0.36 

+0.1 2 44 +015 

-0.1 -0.16 
+0.1 1 fiQ +005 

-0.0 -J^-ua „o.03 

+1-5 1 fifi +2.34 

-1.7 ^-"^ -0.00 
+0.4 o o„ +0.66 

-0.4 -•3-°^ -0.18 
+1.1 1 qy +5.08 

-0.3 -J^-J' -2.63 
+0.1 q 90 +0.08 

-0.1 -J-^u _Q 33 

-4 

+ 0.2 r, Qr, +0.08 

-0.2 ^-"^ -0.01 
+0.4 r, .f, +6.40 

-0.7 -^•'il' -1.60 

+0-1 .3 64 +2-44 
-0.2 '^■"^ -0.36 

+0.1 1 on +0.07 

-0.2 -^-"^ -0.05 
+0.1 r, ooc- +0.164 

-0.4 -o.zoo _o.329 

+ 0.3 r, oqc; +3.16 

-1.1 U.»db _o 00 

+0 564 +0 043 

-0.0 "-Jut -0.082 

+0.1 q a, +0.16 

-1.2 "^-"^ -1.01 

+0.0 9 oq +0.27 
-0.1 -0.20 

+0.2 9 +0.11 

-0.3 ^-^^ -0.16 
+0.1 1 09 +0.05 

-0.1 ''■■"'^ -0.1 
-4 

+0.1 q, 4 +0.33 

-0.5 -•^•J-'i -0.08 
2.06 

+ 0.5 q +0.55 

-0.0 -'J-f'' -0.06 

+0.1 9 +0.20 
-0.2 "° -0.08 

2.2 

+0.4 q oq +3.87 

-0.1 -J-f-J -0.07 
+0.2 q f.~ +0.33 

-0.1 -•J-U' -0.33 



TV TTJ T/'AQQ C 

MrtlvUooO 


22.9 


■NrPPn9fi9 

IN *^\_yUZO.Z 


22 3 




91 A 




21.9 


^MAbXJUloUUzDO 


21.4 




21.9 




23.0 


IN GO 11 94 


22.2 


■\T/^ 1 on 1 

JNCjClzyl 


21.2 


MHKUblo 


20.7 


hjhOoni-CTUlo 


22.3 


100450 


21.9 


UGC03973 


21.1 


MHK07U4 


23.0 


JNGCo227 


21.4 


MESS1ER066 


21.7 


NGC3976 


21.6 


NGC4013 


22.8 


ARP244 


20.9 


NGC4253 


20.9 


NGC4414 


21.1 


NGC4490 


20.7 


NGC4559 


20.7 


NGC4593 


21.3 


MESSIER104 


20.9 


NGC4631 


21.1 


NGC4666 


21.8 


MESSIER063 


23.4 


NGC5194 


23.5 


MESSIER083 


22.1 


IRASF13349-t-2438 


21.1 


NGC5548 


20.8 


NGC5775 


21.6 


MRK0509 


20.2 


IC5169 


21.3 


NGC7314 


20.9 


NGC7469 


20.8 


ESO148-IG002 


21.3 


NGC7771 


21.1 
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Table 8. Parameters of the iron Ka line, be it narrow or broad, if it is present. Column (1) Galaxy name; Column (2) the EW of the 
narrow line at 6.4 keV modelled by either trans or pexmon; Column (3) the energy of the line constrained to be between 6.2 and 6.6 keV, 
not associated with any other component; Column (4) the gaussian width of this component; Column (5) the EW of this component. 



Name EW6.4 EpeKa O-peKa EWpeKa 

(1) (2) (3) (4) (5) 



MRK0335 


31+8 


6.51+°-?? 

"•^-^—0.11 


0.42+°1t 
" — 0.13 


134+1^ 

^ — 38 


NGC0017 


214+2f0 

— 214 








NGC0262 


28+^ 

— 7 














n ,1+0.33 


,40+90 


NGC0424 


576+11^ 

— 187 


6.51+11"'? 


< 1 


109+?24 

^^—109 


NGC0526A 


26+'^ 


6.5+°-l 

— 0.1 


0.34+°f 

— 0.1 


—36 


NGC0513 


116+11 

— 48 


6.2+°-2 

— 0.0 


< 0.35 


48+1°° 

— 39 


2MASXJ01500266 




n 79 + 0.26 
^- ' ^-0.23 


28+"-*^ 

^■^"-0.28 


715+628 

' J^^-448 


NGC1052 


124+^3 


6.51+°?? 

— U.ii 


— u.oo 


23+?? 


NGC1068 


281 


6.6 


0.18 


396 


ARP118 


60+21 


6.47+°i| 

— U.Uo 


0.04+°i6 

— U.U4 


40+« 

— o4 


MCG-02-08-039 


220+11° 








NGC1194 


194ll« 


6.6in 


0.26l°;?® 


23911^1 


NGC1313 


28l28 
-,(,-,+163 




_ 




NGC1^2n 

i 'i V_ J V_y X o ^ u 


6 5+°l 
°-°-o.i 


06+° °^ 


239+209 


NGC1365 


78 








NGC1386 


I7IO+37I 


_ 


_ 


_ 


3C120 


19+^ 


6 43+°-°5 


15+°-22 

^•-^^—0.05 


47+?i 
^' -16 


MRK0618 


30+iJ 

— oU 


6.2+°-2 
— u.u 


< 1.77 


45+1?^ 

— 4i 


NGC1808 


lO+fg 

— iU 








ESO362-G018 


257+®^ 


6.54+°-°5 

— U.Uo 


o.i2i°:?i 


108+if 


2MASXJ05210136 


— 56 






134+^2 


2MASXJ05580206 




6.37+°-l! 

— 0.17 


0.11+°1? 

— U. 11 


45+!^ 

—40 


IC0450 


63+?° 








UGC03973 


6-43i;i;;^ 


o-iitro7 


2051^1 


NGC2655 


140+1!° 

— 140 








MRK0704 


30+23 


6-29+^9 
— u.uy 


0.15+°12 

— U. io 


132+^1 


UGC05101 


127+1^? 


6.55+°°5 

— U.U4 


< 0.29 


170+?59 

— ioo 


NGC2992 


32+^ 


6.33l°12 


0.33l°?^ 


60126 


MESSIER081 


38+11 

— ii 


6.6+Sl 

— U. i 


< 0.13 


23+?4 

— i4 


MESSIER082 




6.48 


< 0.01 


44 


NGC3079 


411 + 173 
— 173 


6.6in 


0.07l°J? 


305l?i8 


NGC3147 


n+_V, 


6-5lS;l 




1621?!^ 


NGC3227 


174+^^ 








NGC3516 




6.41+°-°l 


0-12+'^Jl 


132+19 


NGC3690 


_ 


6.6l°° 


o.3ii°:?i 


913+38° 
'^ — 346 


NGC4051 


117+« 


6-55l°?t 


—0.06 




NGC4151 


174+12 


6.5in' 


< 0.07 




NGC4253 


45+8 






33l^ 


NGC4388 


1441-, 








MESS1ER058 


90111 


6 6+°° 
"•°-o.i 


2+°-2 


206+128 


NGC4593 


108+26 








NGC4631 


492+431 








NGC4666 


ir;r;+100 


6 47+0-04 
°-^'-0.21 


< 0.4 


234+229 
^■^^-170 


UGC08058 


451011^^? 


6 49+0-11 
°-^^-0.13 


24+°-26 
'^■^^-o.ii 


cr,o+967 
J^'J_230 


NGC4968 








MCG-03-34-064 


139lf, 


6 2+°-l 


S6+°-3'' 


254l«l 


NGC5194 


470+160 


6 5+0-1 


< 0.2 


176+pg 


ESO383-G035 


23lg 


6.21E!:; 


51+0-5 


185_2g 


IRASF13349+2438 




6.45l°-i5 


61+° '^'' 


774l2fg 


NGC5256 


36211- 








MRK0273 


202tl°t 






-1 07+95 

l'"-92 


IC4329A 


3611 


6 42+0-08 

"•^^-0.05 


3+°--' 


79 + 16 
'^-13 


UGC08850 


1601^1 


6 6+°° 
"•°-o.i 


06+° " 

U.UD_(j Qg 


1321S05 
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Table 9. Continued. 



Name EW6.4 EpeKa O-peKa EWpoK, 



NGC5506 


-i-Q 


6 


43+ 


„ LA 

0.28Iq 


.33 

.08 


116 


+33 
-29 


NGC5548 


22_y 


6 


,1+0.1 
•4Io.i 


n 1 r+O'. 

0.15+^ 


.18 
.05 


61+28 


2MASXJ151 15979 




6 


9+0.2 
■^-0.0 


r> Qi +0. 

0.31+0 


,46 
.14 


,o-+160 
^•5' -106 


NGC6240 


312+67 










134+60 


NGC6552 


1490+^1° 

— 510 














AM1925-724 


194+239 














MRK0509 


39+^ 


6 


c;a+0.02 
•'^''-0.04 


0.00+J] 


()6 
.01 


17+ 


8 
6 


ESO286-IG019 


491+718 

^^-^-491 














NGC7172 


36tl 


6 


49+0.01 


< 0.01 




361 


10 
8 


NGC7213 




6, 


QQ + 0.02 

"^^-0.02 


0.05l°; 


,03 
.05 




3C445 
















NGC7314 


16t« 


6, 


45+0.12 


0-35l°: 


38 
.11 


121 


+42 
-35 


NGC7469 
















NGC7479 


768tf/9 














NGC7582 


370^29 


6 


g+0.0 
■"-0.1 


0.29^0 


,19 
.18 


269 


+80 
-173 


NGC7674 


070+173 

^'^-173 














NGC7771 




6, 


Qo + 0.22 

•'^^'-o.is 


o.i8i°: 


27 
.18 


487 


+876 
-365 



Table 10. Details of the Compton thick sources found in this sample. Column (1) Galaxy name; Column (2) the directly measured 
ArH(cm~^); Column (3) the EW of the Fe Ka line (eV); Column (4) the reflection fraction determined from the pexmon model. 



Name 


loglo(WH) 


EW 


R 


(1) 


(2) 


(3) 


(4) 


NGC1068 




677 


10.3 


NGC1320 


24.3l,V.3 


940+?-^; 


40 4+^*-^ 
^"■^-39 


NGC1386 


2o.7i;;;i 




145li^? 


NGC1667 








NGC3079 


20.71°-^ 


— 525 




NGC3690 


20.8l°j 


9,0+380 




MESSIER088 


21.31°:^ 




17 6+2^-2 

-^'•°-17.6 


NGC4968 




451011^0 




NGC5194 


20.6+°l 


646lf,f 


CI 0+31.6 


2MASXJ15504152 


21 6+°-* 




1 7+59.5 
' -8.87 


NGC6552 






VI 9+83.2 
-^-^^-47.7 


NGC6890 


21 0+°-^ 




52.ll«.2 


NGC7479 


24 3+°-^ 


768lS^ 


0.29tlf 


ESO148-IG002 






41.6132.6 


NGC7582 


20.6 


639li?I 


107 


NGC7674 




079 + 173 
^'^-173 


,43+120 

^^■^-148 



